Jet-parton inelastic interaction beyond eikonal approximation 
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Most of the models of jet quenching generaUy assumes that a jet always travels in a straight 
eikonal path, which is indeed true for sufficiently hard jet but may not be a good one for moderate 
and low momentum jet. In this article an attempt has been made to relax part of this approximation 
for 2 — 3 processes and found a (15-20)% suppression in the differential cross-section for moderately 
hard jets because of the noneikonal effects. In particular, for the process qq' — >■ qq'g in the centre of 
momentum frame the scattering with an angle wider than ±0.527r is literally forbidden unlike the 
process gg — >■ ggg that allows an angular range ±7r. This may have consequence on the suppression 
of hadronic spectra at low transverse momenta. 



INTRODUCTION 

Constituent quark number scaling of elliptic flow and 
jet quenching are supposed to be the most prominent 
signatures that favour the partonic degrees of freedom 
in the deconfined QCD matter. Overwhelming evidences 
of both signatures coming from dedicated heavy ion ex- 
periments viz., STAR and PHENIX@ RHIC BNL [l|,[3], 
ALICE @LHC CERN estabhshed the fact that, the 
primordial hot-soup of nuclear matter, produced in those 
experiments, indeed contain partonic degrees of freedom, 
before freeze out to hadrons in later stage, instead of 
hadronic degrees of freedom throughout. Observation of 
strong suppression of inclusive yields of high momentum 
hadrons and semi-inclusive rate of azimuthal back-to- 
back high momentum hadron pairs relative to p-p colli- 
sions, are expectations from jet quenching. Both of them 
are extensively explored in collisions of Au-Au nuclei at 
= 200 A GeV in RHIC. ALICE, the dedicated heavy 
Ion collider experiment at CERN, seems to appear as a 
factory of Jets. Evidence for jet quenching has also been 
observed recently in Pb-Pb collision at ALICE Q. 

There are a few well known models in the literature 0- 
S [lll - flij that aim to quantify energy loss (mainly ra- 
diative) and jet quenching phenomena within perturba- 
tivc QCD. Most of these formalisms for energy loss of a 
high momentum parton through gluon radiations, have 
some common technical approximations in order to make 
the calculation simpler. These approximations are imple- 
mented both at the level of single emission kernel calcula- 
tions and at multiple gluon emission estimation schemes. 
Main kinematic approximations at the level of siriqle 
emission kernel, are listed below (also see Refs. hjHia 
for a comprehensive discussion): 

• Eikonal parton trajectory I : Leading parton is hav- 
ing energy E {i.e. = E and p± = 0, by definition 
to start with) is much larger than the transverse 
momentum exchanged gluon with the medium, 
E ^ q±, so that it does not give sufficient trans- 
verse kick to defiect the parton from straight tra- 
jectory along z axis. To relax this approximation 



it is therefore important to keep track of the terms 
of 0{q±/E) in the formalism. 

• Eikonal parton trajectory II : Energy of the lead- 
ing parton is sufficiently high, E ^ kj_ {k± being 
transverse component of the emitted gluon) so that 
it docs not get enough transverse kick also from 
emitted gluon. This ensures that the leading par- 
ton is in eikonal trajectory. However it does not 
fix any definite direction for gluon emission, which 
requires comparison of k± with longitudinal com- 
ponent kz or with energy lu of the emitted gluon. 
Therefore, it is important not to neglect terms of 
0{k±/E) in the formalism to relax this approxi- 
mation in the jet .studies. 

• Soft gluon emission : One often uses the addi- 
tional approximation that the gluon energy is much 
smaller than the leading parton energy uj ^ E. 
When X is the fraction of energy carried out by 
the emitted gluon, i.e., x = w/E, this approxima- 
tion ensures that a; — > 0. When x is typical light 
cone variable, defined by the fraction of light cone 
{+)ve momentum carried out by the emitted gluon, 
i.e., X = k^/p'^ = {k±/ y/s)e'^ , this approximation 
ensures a; — ^ 0, only in the mid rapidity and back- 
ward rapidity regions (— oo > rj > 0) but not in the 
forward regions where rj could be a large positive 
number. 

• Small angle/ collinear gluon emission : Energy of 
the emitted gluon lo is much larger than its trans- 
verse momentum A;^, w 3> fc_L and w ~ fcz. For 
any 2 — >■ 3 process, without loss of generality, one 
can take k_\_ = u) sin 9g and = OJ cos 9g , where 9g 
being the angle between direction of propagation of 
leading parton and direction of emitted gluon. This 



particular approximation therefore implies 



0. 



At this point it is worth mentioning that soft gluon 
emission approximation is a broader approximation as 
it automatically encompasses the eikonal parton trajec- 
tory II approximation, because energy w should always 



be more than the transverse momentum fcj^ for a massless 
emitted gluons. 

In this article we make an effort to relax the eikonal 
parton trajectory I approximation in some extent for the 
radiative/inelastic process 2 — 3. Investigation of all the 
matrix elements in 0{al) for 2 — > 3 radiative processes 
have been done keeping terms up to O {t/s). The first or- 
der in eikonal expansion, i.e., O (t/s), is termed here as 
noneikonal, since the calculation is performed in Feyn- 
man gauge with Mandelstum variables instead of most 
extensively used light cone gauge with light cone vari- 
ablcfQ. We have neglected terms of 0{t^ /s^) and higher 
orders. Throughout our study we also do not use any ap- 
proximation related to small angle/ collinear gluon emis- 
sion. Nevertheless, we have used soft gluon emission ap- 
proximation which automatically include Eikonal parton 
trajectory II assumption. 

Here we relaxed the kinematic constrains associated 
with eikonal propagation. Nevertheless (non)eikonal 
propagation is directly related to space-time tracks that 
the partons take when going through the medium. While 
an eikonal path is a straight line, a non-eikonal path 
is generically any trajectory that have deviation from 
straight line and obviously always somewhat longer. This 
longer path, taken by the parton may matter, in the soft 
sector when the problem is embedded into a hydrody- 
namically evolving density distribution We note 

that there is transport model, although uses the cross 
sections which are kind of eikonal, takes all the trajec- 
tories into account, so particles do not need to move on 
straight lines |2l| . 
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1. Five tree level Feynman diagrams for the process 

qq'g- 




FIG. 2. Non-zero angular deviation from eikonal trajectory. 
Angle between incoming and outgoing momentum of projec- 
tile is 6q and direction of emition of gluon with that of incom- 
ing momentum of projectile parton is 9g. 



with 



s-t-t-fu + s' + t' + u' = 



(2) 



INELASTIC QUARK-QUARK SCATTERING 

(qq' qq'g) 

The process qq' — qq'g (prime denotes different quark 
flavour) in 0{al) appears in five t channel Feynman dia- 
grams, which are shown in the Fig. [T] (see also Fig. [5] for 
other details). Note that ki and /c2 are momenta of the 
different quark flavours in the entrance channel whereas 
^3 and fc4 are those for exit channel and is that of the 
emitted gluon. Scattering angle between ki and kg is 9q 
whereas 9g is the angle between direction of emission of 
soft gluon ks and direction of incoming projectile quark 
ki . We now define the relevant Mandelstum variables for 
this 2 — >• 3 process as 

s = (/^l+/^2)^ s' = (fc3 + fc4)^ 
u = (fci - fc4)^ , u' = {k2 - ks)'^, 

t=(fcl-fc3)2, t' = (fc2-fc4)', (1) 



^ After connecting t to qj_ and s to in the centre of momentum 
frame, term of O (t/s) ensures the relaxation of the approxima- 
tion E ^ q± (eikonal parton trajectory T). 



When the emitted gluon is soft (fcs 0) compare to 
other external legs, one can assume : t' ^ t, s' — > 
s, u' u and we can express the transverse component 
of the momentum of the emitted gluon in the centre of 
momentum frame of fci and k2 as 

2 _ 4(fci • kz)[k2 ■ k^) 
s 

_ js + t + u){s + u' + t') 

= (1±1±^, (3) 

s 

The hierarchy among various scale of momentum, em- 
ployed in the present work is stated as 

^^, £; > VW, q±'>uj>k^>md, (4) 

where is the Debye screening mass acts as an in- 
frared cut-off. We note that the above hierarchy relaxes 
the aproximations ^/s, E 3> 1^ [Eikonal parton tra- 
jectories I ) and w ^ ki^ [small angle/ collinear gluon 
emission ) in some extent. 
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Matrix Elements, Amplitude and Cross-section 

The gauge invariant amplitude summed over all the fi- 
nal states and averaged over initial states for the process, 
qq' qq'g, is 



\M 



qq'~^qq'g\ 



(5) 



where i and j run from to 4. We note that the 
index represents the diagram where the soft gluon 
emits from the exchanged gluon line whereas the indices 
m = 1,2,3,4 represent the diagrams where emission of 
the soft gluons are being from external fermion lines hav- 
ing momenta km (see Fig. [2]). Equation ([5]) contains total 
fifteen terms in which there are five self interfering 'gen- 
uine amplitudes' for i = j and ten cross interfering 'inter- 
ference amplitudes' associated with gluon emissions in- 
volving two diagrams for i ^ j. Here we have extensively 
used the FORM, REDUCE and CalcHep Programs [Ij. 
Results are given below up to 0(1/A;^) and 0{t/s), for 
soft gluon emission. 



Genuine amplitudes 

By genuine amplitudes we are referring those are com- 
ing from each diagram that interfers with itself. In the 
Feynman gauge □ all of them vanish within soft gluon 
emission approximations and in ©(l/fc^): 
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0; M^. 



Ml, = 0; Mlo = 0. (6) 



However, they may contribute in 0(1), 0{k'^), 0{k'^) 
etc, and in 0{t'^/s^) and higher orders. All of them can 
safely be neglected in the soft emission limit as our aim 
is to go beyond the eikonal approximation- 1, E q±. 



Interference amplitudes 

Within the approximations employed above the am- 
plitudes corresponding to the matrix elements of (1 (g) 4) 
and (2 (E) 3) are identical in the leading order (©(l/fe^)) 
as well in 0{t/s) and given as 



Ml 



M 



7128 
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8 27 



1 + 2- 

s 



where 0(t/s) is purely noneikonal (i.e., the first order in 
eikonal approximation) in nature as noted earlier. Also 



the amplitudes for (1 (g) 2) and (3 (E) 4) are identical and 
the contribution is obtained in ©(l/fc^) and 0{t/s) as 



Ml2=Ml 
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t 



Both (1 (E) 3) and (2 E) 4) are also identical within the 
employed hierarchy. However, they do not contribute 
in leading order but only in 0{t/s). Hence, in Feyn- 
man gauge the contribution from the interference be- 
tween initial state and final state radiations, is exclu- 
sively noneikonal in nature and given as 



ML 



1 128 



1 



4 27 ^ t2 k 



1 t 
2l 



Finally any diagram interfering with 0, i.e., (0 E) I) 
with / = 1,2,3,4, docs not contribute in 0{l/k\) but 
contributes in {l/k±\/t) and higher orders. In the limit 
\\/t\ ^ q± ^ w, amplitudes of 0{l/k±\/t) are sublead- 
ing, in comparison to 0{l/k'^). Therefore, all of these 
{MIq, M20, MIq and MIq) do not contribute within 
the approximation used in this work. 

Here we note that color coherence pattern between ini- 
tial and final state radiation in the presence of a QCD 
medium derived in light cone gauge shows angular distri- 
bution of the induced gluon spectrum is broadly modified 
when one includes interference terms [2^ H^l . In Feyn- 
man gauge within our hierarchy only interferance terms 
are contributing to the amplitude. 

The gauge invariant amplitude for the process, qq' — > 
qq'g, can now be obtained by summing all the subampli- 
tudcs as 



--12g^ \M 



qq'-^qq' lekni 



1 

kl 



17t_ 



where the two body amplitude is given as 



^11' \eknl 



8 4S- 

9^ ¥■ 



(8) 



The three-body amplitude in (jT)) for the inelastic pro- 
cess, qq' — > qq'g, contains the two-body amplitude for 
the elastic process, an infrared factor for the emission 
of a soft gluon and a noneikonal correction factor. The 
expression in ([7|) will lead to q± ^ k± limit of the Gu- 
nion and Bertsch formula 25|. If the emitted gluon is 
much softer than others it can then be regulated by the 
Debye screening mass, m^. Terms within the parenthe- 
sis in ([7]) would correspond to noneikonal correction over 
the eikonal Gunion Bertsch formula. Eq. ([7]) is complete 
upto C(l/fc^) and 0{t/s), for emission of a soft gluon in 
the process qq' — >■ qq'g. Similar investigation have been 
done earlier for the process gg — > ggg [iS-Eii. 



Rutherford scattering beyond eikonal approximation 



^ In light cone gauge amplitudes coming from diagrams that in- 
volve gluon emission from target partons can only be neglected, 
others are not. 



It is interesting to note how noneikonality gives way to 
probe beyond Rutherford scattering limits. In the centre 
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eikonal + first order non-eikonal (Eq. 1 1) 




-1 1 

Scattering angle 0q (Radian) 

FIG. 3. Dependence of inverse cross-section cr~^ on scattering 
angle 9q, all prefactor that are independent of 6q taken to be 
unity for convenience. 



of momentum frame for a typical 2 — )■ 2 (or even in case 
of 2 — )• 3 when fifth particle is ultra soft!) process 



- = - snr — 
s 2 



(9) 



The eikonal cross-section (aekni) is directly connected to 
the Rutherford scattering scattering cross-section as 



1 



(10) 



When one relaxes the eikonal approximation then the 
cross-ection can be written as 



tT„eOC 



t2 



1 



nt 



1 



sm 



1 



17 



2 '^9 



(11) 



which puts a restriction on the scattering angle, 9q. In 
Figl3] we have plotted inverse cross-section (ct~^) for 
both eikonal and noncikonal case. Even though both 
behave identically with a similar plateau in the small 
angle scattering but noncikonal cross-section has a very 
sharp fall in comparison with the eikonal one for large 
angle scattering. As seen the noncikonal inverse ma- 
trix element is bounded by the scattering angle, 9q = 
±2 sin~^(3/vT7) ~ ±0.527r, in centre of momentum 
frame, in contrary to that of eikonal one having a natural 
bound of ±7r. This indicates that the back scattering is 
forbidden for the case of qq' — > qq' g when the emitted 
gluon is soft. 

Cross-section in the first order in eikonal (viz. ,noneikonal) 
approximation 

The cross-section for the process qq' qq' g can be ob- 
tained as 



<Jqq'^qq'g 



1 



n 



2 — 3 



\Mq 



I (27r)*J*(fci + k2 



h). 



In the centre of momentum frame, ki + k2 = ka + k4 + ks = 0, and one obtains 

1 f d^ks 1 d^k^ 



'99 -»<?<? 9 



2.S J {2tt)^2E3 {27t)^2E4 {2TTf2E5 



\Mqq.^qq.g\ {2tt)H{Ei + E2 - E3 + E^ + E5) 



1 

2^ 



1 

2^ 



1 1 



2 (27r) 



1 1 



dq\dq^ 
E. 



1 



{2itY2E4, 



1 1 



4 (27r)2 



X 1 Vg/^^ I ^ + y ^ ) (2^)''^(^i +E2-E, + E,+u:) 



2 (27r) 



El 



{2nf2Ei 



1 1 



4 (27r)2 



dk\dri 



17 d 



1 + ^) iT(i-ft)(2-r, 



(12) 



(13) 



where we have used rapidity ?/ = — In [tan (0^/2)], q 
ki — ks, = q sin^q. The cross-section contains fac- 
tors, having term like q\/ s, that are responsible for non- 
eikonal effects. 



In thermal medium taking debye mass as infrared reg- 
ulator, the differential cross-section can be expressed as 



dqj^dkj^drj 



2CACqq>a 



Tab 



ill 



,2\2 



-(14) 
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where Ca = 3 and Cgg' =8/9 are Casimir factors, and 
Tab = CaCb: with various factors arc, cxphcitly, given as 



Ca 
Ct 



9 s 



(15) 



The factor (^a comes from cikonal part of the matrix ele- 
ments, and Cb is the noneikonal factor originated from 
noneikonal part of matrix elements. The differential 
cross-section for the process qq' — > qq' a a.s given in 
Eq. ([T4| correctly reproduces the result of [29[ in the limit 
^ ^1 and in the eikonal limit, y/s,E ^ g^, as all 
the noneikonal factors, viz., (^a and Cb become identically 
unity and so as Tat- 



INELASTIC GLUON-GLUON FUSION (,g,g ^ ggg) 

The three gluon production via gluon-gluon fusion 
gg — > ggg is extremely important in the context of 
heavy-ion phenomenology. For a sequence of events: hot 
glue scenario of glasma field, thermal equilibration, gluon 
chemical equilibration in later time, parton matter vis- 
cosity, radiative energy-loss of high energy partons jet 
propagating through thermalised QGP, this process plays 
a crucial role. Matrix elements for the process gg ggg 
have been computed up to 0{t^/s^) in [27|. Considering 
0{t/s) result it is now straightforward to evaluate the 
differential cross-section for this process in first order in 
eikonal approximation as 



dq'^dk^dr] 



- 2CaC. 



Tab 



1 



{q\ + to2)2 fc2 ^ , 



(16) 



where Cgg ~ 9/2 and the factor. Tab = CaCb, with its 
various components 

„2 ^ -2 



(17) 



The factor comming from eikonal part of matrix el- 
ement (a is same for both processes qq' qq' g and 
gg — > ggg whereas the noneikonal factor C^^ is different. 
This noneikonal factor, obviously, does not put any re- 
striction on the scattering angle [9 = ±7r) for the process 
gg — >■ ggg, and allows it to go in full natural range ±7r as 
compared to the process qq' — > qq'g. 

Unlike gg — )■ ggg where a Park- Taylor type formula 
(30| is available to compute the matrix element, the com- 
putation of matrix elements up to 0{t/s) is quite cum- 
bersome in case of qg — >■ qgg. Also in this article we 
have performed our study on inelastic quark-quark scat- 
tering but with different flavours. In case of same flavour 
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FIG. 4. Typical estimation of noneikonal factors at T = 
SOOAfeV with a = 0.3 for the process qq' — > qq'g- First order 
noneikonal factors (a = (l + 9±/s) i Cft — (l ~ l^Qx/Qs), 
and the full contribution Tat = (aCb- 
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FIG. 5. Typical estimation of noneikonal factors at T = 
SOOMeV with a = 0.3 for the process gg — >■ ggg. First order 
noneikonal factors = (l -f q\/ s) ,(";, = (l — q\/2s) and 
the full contribution Tab = CaCb- 



qq — ?> qqg things would be a more involved one. We leave 
them for future study. 



RESULTS AND DISCUSSION 

For quantitative estimation of the noneikonal effects, 
we have taken the average value of the momentum trans- 
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fer squared which can be obtained [31 as 




2g'T'\n{ElgT). 



(18) 



(q^) has then been embedded in C,a, Cb to have a quah- 
tative estimation of noneikonal effects over eikonal cross- 
sections. 

In FigU and Fig IS] the first order noneikonal factors: 
Ca, Cb and the full contribution Tab for both processes 
qq' qq'g and gg — >■ ggg, respectively, displayed. It 
can be seen that the noneikonal effects are ~ (15 — 20)% 
over eikonal one for moderately hard jets. However, the 
noneikonal effect gradually becomes mild for very high 
energetic jets. Typically for charged hadrons in the mo- 
mentum range of 8 to 50 GeV, i.e. typical parton kine- 
matics of 16 to 150 GeV, way off the scale of Fig|4] and 
FiglSl indicate that non-eikonal effects are largely absent 
above (10 — 15) GeV parton kinematics. In the liter- 
ature attempts have already been made to address the 
noneikonal propagation of partons for collision/clastic 
processes in a monte-carlo approach 3l| by considering 
full 0{a'^) matrix elements for relevant 2 — )■ 2 processes. 
Eikonal propagation approximation was found to be good 
on the 10% level. Present study also revels that for ra- 
diative/inelastic process, Eikonal parton trajectory / ap- 
proximation seems to be within (15 — 20)% level. This 
approximation should be crude only in soft and moderate 
momentum regimes. We also note that splitting kernels 
for partons produced in large virtuality scattering pro- 
cesses that subsequently traverse a region of strongly- 
interacting matter have been investigated early in the 
literature within effective theory formalism [33| . 

There has always been a quest for large angle radia- 
tions. In the present study we do not assume small an- 
gle/collinear emission approximation either in the course 
of calculating matrix elements or in the calculations of 
kinematics. In the present calculation the kinematic re- 
lation kj_ = Lu sin 0g ensures that one can go safely to the 
limit where fc^ ~ w. 

Most of the pQCD based model calculations for de- 
scribing the medium is not able to account the almost 
ideal fluid behaviour, which seems to be a manifesta- 
tion of long range correlations. It also would lead to a 
large elastic contribution to energy loss [s^ for reason- 
able values of coupling which is not supported by the 
data. Sort of single hard scenario 0, [l^ has been dis- 
cussed here at the level of single emission karnel in which 
gluons are induced by a single hard scattering with the 
mediurn, but many widely used quenching models, for in- 
stance [5| , make the opposite assumption of multiple soft 
interactions with the medium leading to induced gluon 
radiation. Such classes of models can probably be ex- 
pected to have somewhat milder non-eikonal corrections 



due to scattering with the medium. 

The kinematic constraints, E ^ lo k±^ q± referred 
in the literature as soft eikonal approximation that ne- 
glects any change in parton trajectory due to multi- 
ple scatterings but assumes a straight line tragectory 
throughout. The diffusion of partons in a hot and dense 
medium can have an unavoidable link beyond the eikonal 
approximation and it is worth to relax eikonal approxi- 
mation. In this work an attempt has been made to re- 
lax part of this approximation for some of the inelastic 
processes and their differential cross-sections in first or- 
der noneikonal approximation have been obtained. Pri- 
mary estimation indicates 15 — 20% reduction in the 
cross section due to first order noneikonal effect for both 
the processes in the soft and intermediate parton ener- 
gies. These cross-sections naturally reproduce eikonally 
approximated results in the eikonal limit for soft emis- 
sion, i.e., \/s,E 3> qi_ and q±_ ^ k^. QGP produced 
at LHC, where large virtuality scattering processes may 
be dominant one, is seems to be 'less opaque to jets 
than predicted' by constrained extrapolations from RHIC 
341. There are however other views also, for instance 
19 , [i^, [3^ , where another set of constrained extrapo- 
lations show considerable variation in the postdictions 
of RHIC-constrained scenarios with LHC data. Here 
this have been taken as a constraint and cause to dis- 
regard class of models which fail to predict/postdict cor- 
rectly the uprising behaviour of nuclear modification fac- 
tor rather than assigning it as generic surprising feature 
of LHC data. 

Our results indicate some reductions in interaction 
strengths of jets due to non-eikonal effects, in soft and 
intermediate sector. In the soft sector when the prob- 
lem is embedded into a hydrodynamically evolving den- 
sity distribution this could lead to non-trivial effects. We 
also show that wide back scattering with scattering angle 
more than ~ ±0.527r is forbidden in case of qq' — >■ qq' g 
when the emitted gluon in soft. This, however, is not the 
case for gg — > ggg. In future study we intend quantitative 
estimation of this noneikonal effect in jet quenching and 
other consequences in heavy-ion collisions phenomenol- 
ogy 
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